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ABSTRACT 

The large scale matter distribution in three different simulations of CDM models 
is investigated and compared with corresponding results of the Zel'dovich theory of 
nonlinear gravitational instability. We show that the basic characteristics of wall-like 
structure elements are well described by this theory, and that they can be expressed 
by the cosmological parameters and a few spectral moments of the perturbation spec- 
trum. Therefore the characteristics of such elements provide reasonable estimates of 
these parameters. We show that the compressed matter is relaxed and gravitationally 
confined, what manifests itself in the existence of walls as (quasi) stationary structure 
elements with life time restricted by their disruption into high density clouds. 

The matter distribution is investigated both in the real and redshift spaces. In 
both cases almost the same particles form the walls, and we estimate differences in 
corresponding wall characteristics. The same methods are applied to several mock 
catalogues of 'galaxies' what allows us to characterize a large scale bias between the 
spatial distribution of dark matter and of simulated 'galaxies'. 
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dark matter - large-scale structure of the Universe 



Over the past decade immense progress was achieved in 
the investigation of the large scale matter distribution. 
Now the galaxy distribution is studied up to the redshift 
z ~ 3 (Steidel at al. 1996). At smaller redshifts the anal- 
ysis of rich galaxy surveys with an effective depth ~ (200 

- 400)/i _1 Mpc, such as the Durham/UKST Galaxy Red- 
shift Survey (Ratcliffe et al. 1996), and the Las Cam- 
panas Redshift Survey (Shectman et al. 1996), have es- 
tablished the existence of wall-like structure elements as a 
typical phenomenon in the visible galaxy distribution in- 
corporating ~ (40 - 50)% of galaxies (Doroshkevich et al. 
1996, hereafter LCRS1; Doroshkevich et al. 1999a, hereafter 
LCRS2; Doroshkevich et al. 2000, hereafter DURS). The 
wall-like structure elements with a typical diameter ~ (30 

- 50) ft -1 Mpc surround low-density regions with a similar 
typical diameter ~ (50 - 70) ft -1 Mpc. Within the wall- like 
structures, the observed galaxy distribution is also inhomo- 
geneous (see, e.g., Fig. 5 of Ramella et al. 1992), and galaxies 
are concentrated in high density clumps and filaments. 

The galaxies occupying low density regions are concen- 
trated within a random network of filaments. Filaments in- 
corporate ~ 50% of galaxies and are clearly seen in many 



redshift surveys (see, e.g., de Lapparent, Geller & Huchra 
1988). These results extend the range of investigated scales 
in the galaxy distribution up to ~ 100ft _1 Mpc. Further 
progress in the study of the observed large scale galaxy 
distribution could be reached with the 2dF redshift survey 
(Colless 1998; Cannon 1998) and the Sloan Digital Sky Sur- 
vey (Loveday & Pier 1998; Maddox 1998). 

The formation and evolution of structure on large scales 
are investigated in numerous simulations (see, e.g., Cole et 
al., 1997, 1998; Jenkins et al. 1998; Governato et al. 1998; 
Miiller et al. 1998; Doroshkevich et al. 1999b, hereafter 
DMRT). These simulations are performed in large boxes (~ 
350 - 500/i -1 Mpc) and reproduce the main properties of 
the observed large scale matter distribution. In particular, 
they confirm formation of large wall-like matter condensa- 
tions due to a nonlinear anisotropic matter compression on 
a typical scale ~ (20 - 30)/i~ Mpc that is about one half of 
the typical wall separation. 

The statistical characteristics of wall formation are de- 
scribed by an approximate theoretical model (Lee & Shan- 
darin 1998; Demiaiiski & Doroshkevich 1999a, b, here- 
after DD99) based on the Zel'dovich nonlinear theory of 
gravitational instability (Zel'dovich 1970, 1978; Shandarin 
& Zel'dovich 1989). This approach relates the structure 
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parameters with the main parameters of the underlying 
cosmological scenario and the initial power spectrum of 
density perturbations. The impact of large scale perturba- 
tions is found to be important throughout all evolutionary 
stages and some statistical characteristics of structure ele- 
ments - filaments and walls formed in the course of nonlin- 
ear evolution - are directly connected with the parameters 
of these perturbations. Another theoretical model of large 
scale structure formation was discussed in Bond, Kofman & 
Pogosyan (1996). 

The simulated large scale matter distribution does not 
exactly reproduce the theoretical expectations due to the 
influence of some essential factors, the most important ones 
are the small scale clustering and relaxation of compressed 
matter, and the large scale matter flow within sheet-like 
structure elements. Thus, compression of matter along one 
of the transversal directions transforms sheet-like elements 
into filaments, while expansion of matter in both transversal 
directions results in the erosion of pancakes. The disruption 
of walls and the small scale clustering of compressed matter 
substantially accelerate the relaxation and are responsible 
for strong matter concentration within walls. This is ap- 
parent from the isotropy of velocity dispersion within walls 
noticed in DMRT. 

The combined influence of these (and other) factors 
complicates the statistical description of the large scale mat- 
ter distribution at late evolutionary stages, what is typical 
for the final evolutionary stages of the standard COBE- 
normalized CDM (SCDM) model with fi m = 1. For low 
density models, such as the open CDM (OCDM) model and 
the ACDM model with 51a > f2 m , the situation is not so 
complex, and some statistical characteristics of structure can 
be successfully compared with the approximate theoretical 
expectations. 

The investigation of wall-like massive structure ele- 
ments is more promising in this respect because walls repre- 
sent the first step in the process of structure formation and, 
so, hold more information about characteristics of the ini- 
tial matter flow. Such walls are observed as superclusters of 
galaxies similar to the Great Wall (de Lapparent et al. 1988) 
and the Pisces-Perseus supercluster (Giovanelli & Haynes 
1993). In simulations such wall- like structure elements are 
also easily identified because of their relatively high overden- 
sity. Samples of such elements were investigated in DMRT 
and LCRS2. The connection between properties of walls and 
the amplitude and the spectrum of initial perturbations was 
discussed in DD99, and some of these results can be com- 
pared with measured properties of simulated wall-like struc- 
ture elements. Examples considered in DD99 and DMRT 
had rather illustrative character, but they seem to be quite 
promising. 

Here we will compare more accurately some of the ex- 
pected and measured characteristics of wall-like matter con- 
densations. We concentrate our attention on the physical 
aspects of the formation and evolution of the large scale mat- 
ter distribution in order to better understand these processes 
and the phenomenon of wall-like matter condensations. Both 
theoretical and numerical estimates are inevitably approx- 
imate, but nevertheless, such comparison allows us to test 
the theoretical conclusions, to reveal and illustrate the in- 
fluence of essential factors mentioned above, and to examine 



the abilities of statistical methods used to describe the large 
scale matter distribution. 

These methods allow us to reveal, in particular, some 
differences in characteristics of the large scale matter dis- 
tribution in the real and redshift spaces. Various aspects of 
this problem were widely discussed during the past decade 
(see, e.g., Kaiser 1987; McGill 1990 a; Davis, Miller & White 
1997; Hamilton 1998; Melott et al. 1998; Hui, Kofman & 
Shandarin 1999, Tadros et al. 1999). Here we show that the 
differences between characteristics of walls in the real and 
redshift spaces depend on the basic cosmological model and 
increase during the cosmic evolution. Characteristics of walls 
in the real and redshift spaces are almost identical for the 
low-density models, but they differ more strongly for the 
SCDM model. 

We do not discuss the application of these methods to 
the observed galaxy catalogues, what is a much harder prob- 
lem, due to the strong influence of selection effects and other 
factors. We will consider this problem in the future. 

This paper is organized as follows: In Sec. 2 the basic 
notations are introduced. In Sec. 3 the statistical characteris- 
tics of wall-like structure elements in the Zel'dovich theory 
are presented. In Sec. 4 we consider the methods used to 
measure the required characteristics of matter distribution. 
Our results are presented in Sees. 5 & 6 where they are also 
compared with the theoretical expectations. Sec. 7 contains 
summary and a short discussion of our main results. Some 
technical details are given in Appendixes A. 



2 STATISTICAL CHARACTERISTICS OF 
LARGE SCALE STRUCTURE 

It is generally recognized that the formation of observed 
large scale structure is driven by the middle part of the 
power spectrum, p(k), with 0.2/iMpc _1 > k > O.OlftMpc -1 
(k is the comoving wave number) , and it is weakly sensitive 
to the small and large scale perturbations. In many publi- 
cations authors use an artificial smoothing of the spectrum 
to describe this process (see, e.g., Bardeen et al. 1986, here- 
after BBKS, Coles et al. 1993). However, as was shown in 
DD99 it is possible to avoid this artificial smoothing if the 
process of structure formation is described in terms of the 
displacement, Si(q), and velocity rather than density field. 

Indeed, in contrast with the density field, the statisti- 
cal characteristics of displacements are weakly sensitive to 
the small and large scale perturbations and are reasonably 
well described by the middle part of the initial power spec- 
trum. Even the strong nonlinear matter clustering does not 
significantly influence the main characteristics of displace- 
ments and, so, such (approximate) description of structure 
holds during long period of structure evolution. Of course, 
this approach cannot describe the formation of gravitation- 
ally confined walls and their disruption into a system of high 
density clouds. 

Bearing in mind these comments we will describe 
the structure parameters using characteristics directly con- 
nected with the displacement. One of them is the large scale 
amplitude of perturbations measured by the dispersion of 
displacements, 

1 

CTs2(z)= 2^ J P( z >^ dk > ( 2A ) 
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Other convenient parameter is the coherent length of the 
displacement and velocity fields, l v , expressed through the 
moment m_2 of the initial power spectrum, p(k). Suitably 
defined coherent length l v provides simple expressions for 
the correlation functions of these fields and the basic char- 
acteristics of the large scale structure (DD99 and Sec. 3). 



2.1 The Zel'dovich approximation 

The Zel'dovich theory connects the Eulerian, n, and the 
Lagrangian, qi, coordinates of fluid elements (particles) by 
the expression 



r l = (l + ^)- 1 [g l -B( 2 )S' l (q)], 
S i (q)=9$(q)/a gi , 



(2.2) 



where z denotes the redshift, B(z) describes growth of per- 
turbations in the linear theory, and the random vector Si or 
the random potential <E> characterize the spatial distribution 
of perturbations. The Lagrangian coordinates of a particle, 
qi, are its unperturbed comoving coordinates. 

The velocity of a particle can be found from (2.2) as 



dri 
~dt 



H{z) 
l + z 



[qi-(l+P)B(z)Si(q)}, 



l + zdB(z) 



B 



dz 



(2.3) 



H(z) = H ^n m (i + zf + (i - n m - n A )(i + z) 2 + n A , 

where H is the Hubble constant (Ho = 100 • h km/s/Mpc). 
Analytical fits for the functions B(z) and (5(z) were given in 
DD99. Approximately, at z <C 1, we have 



5(0) = 1, 0(0) 



2.3Q. 



1 + 1.30, 



(2.4) 



2.2 Main structure characteristics for the 
CDM-like power spectrum 

The standard CDM-like power spectrum with a Harrison 
Zel'dovich large scale asymptote 



p cdm (k) = A(z)kT 2 (k/k ), 



r/iMpc~ 



(2.5) 



l-7p 7 

prel 



can be taken as a reasonable approximation of the initial 
power spectrum used in Zel'dovich' theory. Here A(z) is the 
amplitude of perturbations, T(x) is a transfer function and 
p 1 & prei are the densities of CMB photons and relativistic 
particles (photons, neutrinos etc.). For this spectrum the 
parameters l v and a s are expressed through the spectral 
moments, rrij, as follows: 



kT 2 (k/ko)dk = m,-^k%, 



(2.6) 



i r 

Wo PC 



l (k)dk 



A (Z) ,2 



A(z) 
2-kHI 



f 

Jo 



x A+3 T z (x)dx, m_ 2 



xT (x)dx 



For the CDM transfer function (BBKS) m_ 2 = 0.023, and 
the expressions for the scale l v and the characteristic masses 
of DM and baryonic components associated with the scale 
l v can be written more explicitly as 



6.6 /0.023 



m_2 



h x Mpc, 



(2.7) 



M v = — < p>l% 



2 ■ 10 14 M G 



Here £lt is the dimensionless mean density of the baryonic 
component. The same characteristic scale, l v , as given by 
(2.7) can be used for the structure description as long as the 
Zel'dovich theory can be applied. 

More details can be found in DD99. The same approach 
can be used for other power spectra as well. 



2.3 The amplitude of large scale perturbations 

The large scale amplitude of perturbation as measured by 
A(z) in (2.5) and <j s (2.1) can be successfully used to 
describe the structure evolution in the framework of the 
Zel'dovich theory. As was shown in DD99 it is convenient 
to use - together with a s - an effective dimensionless 'time', 
t(z, Q m , h), 



r(z) 



(2.8) 



which is proportional to the large scale amplitude of per- 
turbations and describes suitably the evolutionary stage 
reached in the model. This 'time' is similar to that used 
in the adhesion model (Shandarin & Zel'dovich 1989). 

As was noticed in DD99 the structure evolution shows 
strong features of self-similarity and is described by uni- 
versal expressions depending on the dimensionless variables 
q/7„ and r. This is a direct consequence of the Zel'dovich 
approximation. 

The 'time' r can be measured by different methods, 
some of which are discussed below. It is sensitive to the sam- 
ple under investigation and to the method of measurement. 
It can be used to quantify bias between spatial distributions 
of different objects, such as, for example, large scale bias 
between distributions of galaxies and the DM component. 

The quadrupole component of the CMB anisotropy, Tq, 
the variance of density in a sphere with radius 8fe _1 Mpc, 
erg, and the velocity dispersion, a ve i are the more often used 
characteristics of the large scale amplitude. All these charac- 
teristics are proportional to each other, but their dependence 
on Q, m and h is different, and they are sensitive to matter 
distribution in different scales. Thus, the quadrupole com- 
ponent of CMB anisotropy characterizes the perturbations 
on scales comparable with the horizon, while the values a ve i 
and as are more sensitive to the matter distribution in mod- 
erate and small scales. 

The connection of these characteristics with o s and r 
can be summarized as follows: 

(i) Using the fits for the CMB anisotropy proposed by 
Bunn & White (1997) we obtain for the flat ACDM and 
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open OCDM models 



2.73 h 2 ^ 2 I m - 2 Tq 



0.023 20^ 



, O.A = 1 — fin 



(2.9) 



^ n 7 o 7,201.65-0.19 lnn r , 



m_2 Tq 
0.023 20/^ 



, Ov = 0, 



where tt denotes the amplitude of large scale perturbations, 
r, measured by the CMB anisotropy. These estimates de- 
pend on the spectral moment m_2 only which is very stable 
and does not change during the considered period of evo- 
lution. But the estimates should be corrected if a possible 
contribution of gravitational waves is taken into account. 

(ii) The amplitude of perturbations, o~ s and r, can be di- 
rectly expressed through the two point autocorrelation func- 
tion as follows: 



lim / dx 1 x^(x), 



(2.10) 



and for the autocorrelation function £(r) approximated by 
the power law 



£(?) = ( r «/r) 7 , r<r 5 , 



we have 



tin) ■■ 



(2-7)(3-7)' 



(2.11) 



(2.12) 



Here is the first zero-point of the autocorrelation function 
and t<: denotes the amplitude r measured by this function. 
The parameter r% is usually found with small precision, but 
for 7 ^1.5 - 1.7, 1 — 7/2 ^0.25 - 0.15 even some variations 
of r£ do not change significantly the final estimate of r. 

(iii) The parameter erg can be also expressed through the 
two point autocorrelation function, £(r), (Peebles 1993), and 
for £(r) approximated by a power law (2.11) we have 



2 

08 



72 



(3-7)(4- 7 )(6-7) V Wh-iMpc 



ro 



r 8 = agF 



(4- 7 )(6-7) 



36.75(2 -7) Kmh^Mpc 



(2.13) 



2- 7 



(2.14) 



(iv) The dispersion of the peculiar velocity of particles at 
small redshifts, z~0, can be written as in the linear theory 
(DD99) 



p- Q m 1535km/s 

<?vei = mov3t, u = IvHofi ~ — - - - ^ 

and for r we obtain the independent estimate 

0~vel 



T V el 



V3 



M0 



(2.15) 



(2.16) 



Here T ve i denotes the amplitude r measured by the veloc- 
ity dispersion. a ve i takes into account also the high veloc- 
ities generated by the gravitational compression of matter 
(in particular, within clusters of galaxies) and, so, it gives 
actually an upper limit of the amplitude. 



3 STATISTICAL CHARACTERISTICS OF 
WALLS IN THE ZEL'DOVICH THEORY 

In both observed and simulated catalogues, at small red- 
shifts, the wall-like structure elements accumulate ~ 50% 
of galaxies and form the skeleton of large scale structure. 
So, investigation of the characteristics of these elements is 
important in itself. It allows us also to obtain information 
about processes of nonlinear structure evolution. In particu- 
lar, we can find two independent measures of the large scale 
amplitude, r. As walls represent the first step of the large 
scale nonlinear matter compression their characteristics can 
be compared with predictions of the Zel'dovich theory. 

In this Section we will consider five characteristics of 
walls, namely, the surface density of walls, m w , defined as 
the amount of matter per unit of wall surface, for example, 
per /i _2 Mpc 2 , the thickness of walls, h w , the wall sepa- 
ration, D sep , the velocity dispersion of matter compressed 
within walls, w w , and the dispersion of wall velocities, a v . 
All these characteristics can be derived from the Zel'dovich 
theory (DD99) and can be found for simulated point distri- 
butions as well. 



3.1 Formation of walls 

Following DD99, we will consider the intersection of two 
fluid particles with Lagrangian coordinates qi and q2 as 
the formation of a wall (Zel'dovich pancake) with the sur- 
face density m w = (n p )\qi — where (n p ) is the mean 
particles density in the sample. In Zel'dovich theory statis- 
tical characteristics of such walls are described by the ini- 
tial power spectrum (2.5) and can be expressed through the 
characteristic scale, /„, the surface density of wall, m w , or 
dimensionless surface density, Qui — /TZ iu /l v /(n p ), and the 
'time', r, introduced in Sees. 2. To do this, the structure 
functions of the initial power spectrum can be used. For the 
standard SDM-like power spectrum (2.5) with the BBKS 
transfer function these functions were introduced in DD99. 

Naturally, the theoretical considerations describe the 
idealized model of structure evolution. Thus, it uses the rigid 
wall boundary though in reality such boundaries are always 
blured. Other important factor is the compression and ex- 
pansion of pancakes in transversal directions. These motions 
transform pancakes into filaments and/or lead to the dissi- 
pation of poor pancakes. They are not so important for rich 
walls but can change the wall surface density by a factor of 
1.3 - 1.5. The small scale clustering and relaxation of mat- 
ter distorts also the measured characteristics of walls with 
respect to theoretical expectations. 

These factors distort the actual power spectrum with re- 
spect to the used one and introduce differences between the 
expected and actually measured parameters of walls which 
however cannot be evaluated a priori. The actual power and 
limitations of this approach must be tested at first with N- 
body simulations. 



3.2 Wall properties in the real space 

3.2.1 Surface density of walls 

The most fundamental characteristic of walls is the surface 
density, m w . The approximate expression for the probability 
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distribution function (PDF) of the pancakes surface density, 
m w , defined as above has been obtained in DD99 in the same 
manner as the well known Press-Schechter mass function. It 
characterizes the process of one dimensional matter com- 
pression and formation of wall-like pancakes as described by 
the Zel'dovich theory. 

For Gaussian initial perturbations and the standard 
CDM-like power spectrum with the BBKS transfer function, 
it can be written as follows: 



N„ 



1 



27TT m \fOw 



q w = 



j{n p ) 



fc eXP (-||) erf (yi|) 



-Q2i r 

'« ' Jo 



(3.1) 



qi- 



N m (q w ) dq w = 1, 



(qw) 



q w N m {q w ) dq w = 8(0.5 + l/n)r, 



— 6.55t 2 



(q w ) = / qiN m (q w ) dq w = 128(0.375 + 1/tt)t* » 887r* , 
Jo 

where (n p ) is the mean particle density in the sample, l v is 
defined by (2.7), r m characterizes the amplitude of pertur- 
bations and the evolution stage of structures, r, as measured 
by the surface density of walls, and qi & q2 are Lagrangian 
coordinates of wall boundaries. This relation was corrected 
for the merging of neighboring walls and this process is de- 
scribed by the erf - function in (3.1). 

These expressions connect r m with the mean surface 
density of walls and allow us to estimate r m from measure- 
ments of (q w ). For other models and/or other distributions 
of initial perturbations the PDFs similar to (3.1) could be 
obtained using the technique described in DD99. 

3.2.2 The wall separation 

We have not been able to find a simple theoretical descrip- 
tion of the wall separation. Nonetheless, taking into account 
the mainly one dimensional character of wall formation, we 
can roughly link the mean measured wall separation, (D 3ep ), 
to the mean surface density of walls, (q w ). 

Indeed, the matter conservation law along the direction 
of wall compression can be approximately written as follows: 

(m w ) w f w (n p )(D sep ), (q w ) ^ f w (D sep ) /l v , 

where f w is the matter fraction assigned to walls. It implies 
that on average a fraction /„, of particles situated at the 
distance ±0.5D aep from the center of wall will be collected 
by the wall. For simulations when the mean wall separation 
is comparable to the box size, L box , we will use the more 
accurate relation 



fdq I Dsep 
l v \ 1 + Dsep/ Lbo 



(3.2) 



The averaging can be performed analytically assuming the 
exponential distribution function for the wall separation. 

The factor fd q defined by Eq. (3.2) characterizes the 
matter fraction assigned to walls as it is determined by com- 
parison of independently measured characteristics (q w ) and 
(D sep ). In turn, difference between f dq and f w characterizes 
the robustness and degree of self consistency of the model 



and the measurements. These estimates are only approxi- 
mations because the wall formation is actually a three di- 
mensional process. 

3.2.3 Velocity of structure elements 

For the pancakes defined in Sec. 3.1 the ID velocity of walls, 
v w , can be found from relations (2.2) and (2.3) as follows: 



1 



qi - q2 



n[u — H(z)r] dq 



(3.3) 



where n is a unit vector normal to the wall. The small scale 
clustering and relaxation of compressed matter does not in- 
fluence the velocities of walls and, so, they are the most 
stable characteristics of the evolutionary stage reached. As 
was shown in DD99 the mean velocity of walls, (v w ), is ex- 
pected to be negligible as compared with its dispersions, a v , 
and the expected PDF of this velocity, N v , is Gaussian for 
Gaussian initial perturbations. 

For the standard spectrum (2.5) with the BBKS trans- 
fer function, and for }„, < 1, the velocity dispersion is re- 
lated to the amplitude of initial perturbations as follows: 

UOT, T v = — , (3.4) 



Uo 



what is similar to (2.16) and also is identical to expecta- 
tions of the linear theory. Here t v denotes the amplitude r 
as measured by the dispersion of wall velocity and uo was 
introduced by (2.15). 



3.2.4 Velocity dispersion of matter compressed within 
walls 

The variance of velocity of matter accumulated by walls, 



qi 



1 r 2 



(z)nr — v w ] 2 dq, (3.5) 



can be found in the framework of the Zel'dovich theory using 
the structure functions described in DD99. As is shown in 
Appendix A, it can be written as follows: 

~ 2 r 2 (l+/3) 2 



■a 1 n 2 qw , 

w wz (q w ,T) w u Q ( — + 



3B 2 



q w < 1, (3.6) 



where /3, u , q w & v w were introduced by (2.3), (2.4), (2.15), 
(3.1) & (3.3). In fact, this function characterizes the mean 
kinetic energy of particles compressed into a wall of a given 
size q w . After averaging over a sample of walls with the PDF 
N m (3.1), in the Zel'dovich theory, we obtain 

wl(r) = r)> « ulr' (7A+ 2 - 2{1 + f3)2 ) . (3.7) 



P 2 

The comparison of the expected mean kinetic energy of the 
compressed particles with the kinetic energy measured in 
simulations characterizes the mean degree of relaxation of 
compressed matter at a given r. 

For richer walls, with q w (q w ), the relation (3.6) is 
transformed into 



12 



Qw, 



(3.8) 



and for such walls, the PDF is similar to (3.1). For a rich 
sample of walls, this relation can be also used for the direct 
measurement of the amplitude r (DMRT; DD99). 



© 0000 RAS, MNRAS 000, 000-000 



6 Demianski et al. 



3.2.5 Wall thickness 

The methods discussed in DD99 allow us also, in the frame- 
work of the Zel'dovich theory, to obtain the expected thick- 
ness of walls along the direction of maximal compression, 
h w . It can be characterized by the thickness of a homoge- 
neous slice with the same surface density. The corresponding 
expression (Appendix A) is 



h wz (q w ,r) « 2I„r^ (1 + z)~ 



(3.9) 



This relation shows that the wall thickness is strongly corre- 
lated with its surface density. After averaging with the PDF 
(3.1) we obtain for the mean thickness of walls 



{h wz ) ^8n- 1/2 l v r\l + z) 



(3.10) 



The degree of matter compression in the Zel'dovich the- 
ory, 5 z (q, t), is characterized by the ratio 



8 Z = 



>qw 



2t 



After averaging with the PDF (3.1) we have for the mean 
degree of matter compression 



2r 



2 

7^ 



1.13. 



(3.11) 



So, in the Zel'dovich theory the averaged degree of matter 
compression is small. 



3.3 Wall properties in the redshift space 

In observed catalogues only the redshift position of galaxies 
along the line-of-sight is known, and therefore the parame- 
ters of observed structures with respect to those found above 
can differ due to the influence of the velocity field. The sta- 
tistical characteristics of walls in redshift space predicted by 
the Zel'dovich theory can be found with the methods de- 
scribed above. This information is not so rich as in the real 
space because in the redshift space, positions of particles 
are determined by their velocities, and, for example, such a 
useful characteristic as the wall velocity cannot be found. 



3. 3. 1 Surface density of walls 

In the real space (Sec. 3.1) the pancake formation was de- 
fined as an intersection of particles with coordinates qi & q2. 
In the redshift space the velocity (2.3) along the line-of-sight 
must be used instead of the coordinate. In Zel'dovich theory 
the velocity dispersion exceeds the dispersion of displace- 
ment by a factor of (1+/3). Hence, this substitution increases 
the wall surface density in the redshift space in respect to 
that in the real space, and now we must use 



Trd 



= frdT = r v /(l + /3) 2 cos0 2 + 



sin ( 



(3.12) 



instead of t. Here the factor f r d >1 describes the more effec- 
tive matter compression in the redshift space predicted by 
the Zel'dovich theory, f3 was introduced in (2.3), (2.4), and (j> 
is a random angle between the direction of wall compression 
and the line-of-sight (0 < <j> < tt/2). Evidently, r rd — r for 
P = 0, so f rd (fi = 0) =1. 

The PDF of wall surface densities in the redshift space 



is identical to (3.1) with a substitution of r r d for r, and now 
for the mean surface density of walls we have 



(q w ) = 8(0.5 + l/n){f? d )T 2 



6.55(/ 2 d }r 2 , 



(3.13) 



l<{frd) = -[2 + (/3 + l) 2 ] < 3.667, 

where r characterizes the evolutionary stage as before. Prob- 
ably, these relations can be used for the description of poorer 
pancakes and earlier evolutionary stages when the influence 
of other factors is less important. 

At small redshifts we must take into account the influ- 
ence of the high velocity dispersion of compressed matter 
generated by the small scale matter clustering and relax- 
ation. The influence of this factor, well known as the 'fin- 
ger of God' effect, is opposite to that discussed above. It 
changes the observed particle position within walls along the 
line-of-sight what blurs the wall boundary and increases the 
thickness of observed walls. It artificially removes the high 
velocity particles from the selected wall and effectively de- 
creases the surface density of walls selected in redshift space 
with respect to the estimates (3.13). 

The impact of this factor can be approximately de- 
scribed by a modification of PDF of wall surface density, 



exp 



erf 



(3.14) 



q w 



8/ 2 



ex p (~8t^) ' w (q™> T > s t. 



and a new normalization of distribution . 

The second term in the square brackets describes the ar- 
tificial rejection of high velocity particles from the wall with 
a surface density q w bounded by a threshold density 8 thr- 
ill this term the exponent gives the fraction of matter accu- 
mulated by the wall in real space for some q w & r, whereas 
W(q w ,T, 5 t hr) is the fraction of high velocity particles which 
are removed from the wall in the redshift space. The func- 
tion W(q w , t, 5 t hr) cannot be found in the Zel'dovich theory 
as it depends on the distributions of particles positions and 
velocities arising due to the small scale clustering and relax- 
ation of matter compressed within walls. 

An other factor which can suppress the expected dif- 
ference of wall characteristics, measured in the real and 
redshift spaces at small redshifts, is the strong matter con- 
densation within structure elements with various richnesses. 
The strong matter rearrangement transforms the continuous 
matter infall on walls into a discontinuous one, increases the 
separation of infalling structure elements, even in the red- 
shift space and, so, at least partly, prevents the erosion of 
wall boundaries. 

These comments show that in the redshift space the 
Zel'dovich theory with the factor / r£ j given by (3.12) 
& (3.13) overestimates the matter concentration within 
walls. Therefore, instead of the factor f r d in (3.13) a factor 
Krd(T,T,l t hr) should be used and the more realistic relation 



6.55k 2 / 



1 K, r d ^ fr 



(3.15) 



connects the amplitude r m with the wall richness q w in the 
redshift space. 
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The actual value of K r d depends on the parameters of 
the cosmological model and on the method of identification 
of walls. The analysis performed below shows that for the 
walls selected in 3D space as described in Sec. 6.1, no growth 
of q w was found, and the parameters q w & r m are connected 
by the relation (3.1) as in the real space. 



3.3.2 Wall separation 

The separations of richer walls is not sensitive to relatively 
small shifts of particle positions introduced by the random 
velocities, but these shifts can result in an artificial merging 
of poorer walls. The influence of this factor can be tested 
with the relation (3.2) as before. 



3.3.3 Velocity dispersion of matter compressed within 
walls and the wall thickness 

In the redshift space the expression for the velocity disper- 
sion of matter compressed within walls in the Zel'dovich 
theory is identical to (3.6) with a substitution of r rd = TK rd 
instead of r, but now it characterizes also the observed thick- 
ness of walls. For walls selected from the 3D sample of par- 
ticles, as is described in Sec. 6.1, we have 



= Vl2w w H 1 , 



"' = "" xl u~ — + 



(3.16) 



(3.17) 



This value exceeds the corresponding real thickness of walls 
given by (3.9). The expected overdensity of compressed mat- 
ter is given by 



0~rd — lv (c[w j hw) . 



(3.18) 



4 MEASURED CHARACTERISTICS OF 

LARGE SCALE MATTER DISTRIBUTION. 

4.1 Core-sampling approach 

The core-sampling approach was proposed by Buryak et al. 
(1994) for the analysis of the galaxy distribution in deep pen- 
cil beam redshift surveys. In the original form it allows to 
obtain the mean free-path between the filaments and walls. 
It was improved and described in detail in LCRS1 where 
some characteristics of the large scale galaxy distribution 
were found for the Las Campanas Redshift Survey. For sim- 
ulated matter distributions as considered here these charac- 
teristics were discussed in DMRT. 

The potential of the core-sampling approach is not ex- 
hausted by these applications, and it could be used to mea- 
sure parameters of the large scale matter distribution dis- 
cussed in the previous Sections. Here we will use this ap- 
proach to obtain the characteristics of the wall-like structure 
component. 

The core-sampling method deals with a sample of points 
(galaxies) lying within relatively narrow cores - rectangular 
and/or cylindrical in simulations, and conical in observations 
- and it studies the point distribution along these cores. For 
some applications the transversal coordinates of points can 
be used as well. To take into account the selection effects, 



which are important for observed catalogues, appropriate 
corrections can be incorporated. The sampling core is char- 
acterized by the size, D core , that is the side of a rectangular 
core or the angular diameter of a conical core. 



4.2 Measured characteristics of walls 

Here we will apply the core sampling technique to the sam- 
ple of wall-like structure elements selected by a 3D-cluster 
analysis (DMRT, Sec. 6.1). This means, the sampling cores 
contain only the particles assigned to walls. Further on, all 
particles are projected onto the core axes and are collected 
into a set of clusters with a linking length Z«nfc- Clusters 
with richness larger than a threshold richness, N m in, arc 
identified with walls within the sampling core. 

The measured wall parameters are sensitive to the in- 
fluence of small scale clustering of matter within walls. For 
strongly disrupted walls and a narrow core, the results de- 
pend on the random position of high density clumps, what 
strongly increases the scatter of measured wall properties. 
The influence of this factor is partly suppressed for larger 
sizes of the sampling core, D core . 

However, the random intersection of the core with a 
wall boundary generates artificially poor clusters. The num- 
ber of such intersections increases proportionally to D core 
what restricts the maximal D core . To suppress the influence 
of this factor a threshold richness of cluster, N m i n , was used. 
If however N mirL becomes too large, the statistical estimates 
become unreliable. For large D core the overlapping of pro- 
jections of neighboring walls becomes also important what 
distorts the measured wall characteristics. 

It is also important to choose an optimal linking length, 
kink, because for small lunk, only the high density part of 
walls is measured, whereas for larger ln n k, again the impact 
of the random overlapping of wall projections becomes im- 
portant. 

The influence of these factors cannot be eliminated 
completely, and our final estimates of properties of walls 
are always distorted to some degree. These distortions 
can be minimized for an optimal range of parameters 
D C ore, N mirl & kink- Practically, these factors do not distort 
the velocity dispersion of walls, a v , which therefore provides 
the best characteristic of the actual evolutionary stage of the 
wall formation. On the other hand, the comparison of results 
obtained for different lunk and D core allows to characterize 
the inner structure of walls. 



4-2.1 Measurement and correction of wall parameters 

The richness of clusters in the core measures the surface 
density of walls, 



N„ 



D 2 

J^co 



(4.1) 



where N m is the number of particles in a cluster. The ve- 
locity of walls, v a im, the velocity dispersion of particles ac- 
cumulated within walls, w a i m , and the proper sizes of walls, 
h sim , are found as follows: 

r w = — y^ri, Vsi m = — y^jui - Hn), 

i=l i=l 
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^ N m 

"L = N _ 1 y^("i - ^ ~ « s im) 2 , (4.2) 
m i=l 

Here r», r w and «i are the coordinates of a particle, of a 
wall, and the velocity of a particle along the sampling core, 
respectively. The wall separation, D sim , is measured by the 
distance between neighboring clusters. 

The parameters m 3im , v 3im , w 3im and h sim as given by 
(4.1) and (4.2) are found along the sampling core and, so, 
are not identical to the parameters discussed in Sec. 3. These 
parameters must be corrected for the random orientation 
of walls with respect to the sampling core. The impact of 
this factor increases the measured surface density, and the 
corrected wall surface density, m c , is connected with the 
measured one by 

m c = m sim cos <j>, < <j> < n/2, 

(m c ) = 0.5(m sim ), (4.3) 

where is a random polar angle between the core and the 
vector orthogonal to the surface of the wall, and the averag- 
ing is performed in a spherical coordinate system. Corrected 
values of the wall velocity and the walls thickness are as fol- 
lows: 

v c = v alm V3, h c = h slm /V3. (4.4) 

In the redshift space the wall thickness is connected with 
the velocity dispersion by Eq. (3.16). The velocity dispersion 
within walls was found to be almost isotropic (DMRT), and, 
so, we will use the measured the actual velocity 

dispersion across walls. 

The measured PDF of the wall surface density, N m (m c ), 
and the mean wall surface density, (m c ), are distorted due 
to the small statistics of rich walls and rejection of poor 
walls with a richness N m < N m in- The correction for these 
distortions can be estimated by comparing the simulated 
PDF with the expected PDF (3.1). 

To do this we will fit the measured PDF to the function 

N m = -^=e~ Xm eri(s/x^), x m = &m?7W " . (4.5) 

The parameter b m describes deviations of measured and ex- 
pected mean surface density of walls (m a i m ), and a m is a 
normalization factor. If the measured PDF is well fitted to 
the function (4.5) then the value 

mt = (m c )/b m (4.6) 

can be taken as a measure of the 'true' mean surface density 
of walls. 

Finally, the mean dimensionless surface density of walls, 
(q w ) and the amplitudes of perturbations, r m & t v , mea- 
sured by the surface density and velocity of wall-like struc- 
ture elements, can be estimated as follows: 

, . _ (m aim ) /M _V^J (47) 

{qw) ~ 2b m l v (n p ) ' Tm -y 6.55 ' Tv ~ no ' ^ 

The small statistics of rich and poor walls distorts also 



the measured wall separation, D aep . The expected distribu- 
tion of wall separations is exponential, and therefore it is 
possible to correct the mean separation using the fit of the 
measured PDF, N sep (D s i m ), to the function 

Nsep — a S e P exp{~bsepDsim/ (Dsim}). (4.8) 

As before, the parameter b aep describes deviations of the 
measured and expected mean separation of walls, and a aep 
is a normalization factor. If the measured PDF is well fitted 
to the function (4.8) then the value 

(D S e P } = (Dsim} /b S e P > (4-9) 

can be taken as a measure of the 'true' mean separation of 
walls. 

5 GENERAL CHARACTERISTICS OF THE 
SIMULATED MATTER DISTRIBUTION 

5.1 Basic simulations 

The theoretical model discussed above describes the evo- 
lution of the DM distribution and, so, should be tested 
with the simulated DM distribution as well. Here we use 
three simulations as a basis for our analysis - the COBE 
normalized standard CDM model (SCDM), a ACDM with 
Qa > fim, and an open CDM (OCDM) model. These models 
were described and investigated with 3D cluster analysis and 
Minimal Spanning Tree technique in DMRT. It was found 
that the ACDM and OCDM models successfully reproduce 
the main observed characteristics of large scale matter dis- 
tribution while the SCDM model demonstrates strong sig- 
natures of overevolution. Here we study these three models 
bearing in mind that only the ACDM and OCDM models 
can be considered as realistic models of the observed large 
scale matter distribution. The SCDM model represents the 
matter distribution typical for a late evolutionary stage. 

The simulations were performed with a PM code 
in a box of (500ft _1 Mpc) 3 with (300) 3 particles for the 
Harrison-Zel'dovich primordial power spectrum and the 
BBKS transfer function. The force and mass resolutions are 
~ 0.9ft _1 Mpc and ~ 1O 11 M , respectively. The point dis- 
tribution in redshift space was produced by adding an ap- 
parent shift to one coordinate due to the peculiar velocity 
of particles. 

Four mock catalogues were prepared on the basis of the 
OCDM model with various degrees of large scale bias be- 
tween the spatial DM distribution and the 'galaxies'. These 
mock catalogues were constructed by identifying randomly 
'galaxies' with DM particles, but with a probability depend- 
ing on the environmental density, thereby identifying more 
particles as 'galaxies' in high density regions (walls). These 
catalogues were investigated also both in real and redshift 
spaces. 

The main characteristics of the simulations are listed in 
Table 1. A more detailed description can be found in DMRT. 

5.2 Large scale amplitude of perturbations 

The evolutionary stages reached in the models under dis- 
cussion can be suitably characterized using the methods de- 
scribed in Sec. 2. The value tt listed in Table 1 characterizes 
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Table 1. Parameters of simulated DM and mock catalogues 



en m nl 


o 


h 


T T 


&S 


r 


7 


T £ 


° vel I V *J 


'vel 












7i _1 Mpc 






km/s 




SCDM 


1 


0.5 


0.68 


1.37 


6.5 


1.9 


0.94 


670 


0.51 


ACDM 


0.35 


0.7 


0.37 


1.11 


6.0 


1.8 


0.34 


554 


0.37 


OCDM 


0.5 


0.6 


0.29 


0.74 


5.0 


1.3 


0.25 


346 


0.23 


mocki 


0.5 


0.6 


0.29 


0.95 


6.0 


1.4 


0.28 


370 


0.25 


mock2 


0.5 


0.6 


0.29 


0.95 


6.0 


1.4 


0.28 


370 


0.24 


rnock^, 


0.5 


0.6 


0.29 


1.24 


7.0 


1.5 


0.33 


374 


0.24 


mocki 


0.5 


0.6 


0.29 


1.61 


8.0 


1.6 


0.39 


404 


0.26 



ro & 7 are the correlation lengths and the slope of correlation function (2.11) in redshift space; 
<TSi T"Ti an d T £ arc the amplitudes of perturbations as given by (2.9) & (2.12), cr ve i & r ve i 
are the velocity dispersion of all particles and amplitudes of perturbations measured by cr vf ,i 
as given by (2.16); 



the large scale amplitude used for the normalization of simu- 
lated perturbations. Other measures of the amplitude, such 
as erg, & T ve i, are sensitive to both the actually realized 
sample of random perturbations (cosmic variance) and to 
the nonlinear distortions of power spectrum produced dur- 
ing the evolution. For the considered mock catalogues these 
measures are also sensitive to the large scale bias between 
the spatial DM and 'galaxies' distributions what allows us 
to characterize it quantitatively. 

The spatial matter distribution and the bias between 
spatial distributions of DM component and 'galaxies' can 
be characterized by the correlation length, ro, and the slop 
of the correlation function, 7, introduced in (2.11). These 
parameters are listed in Table 1 for all samples. Using rela- 
tions (2.12) and (2.13), these values allows to calculate as 
and T£, which are also listed in Table 1. 

The characteristics of correlation function, ro and 7, are 
sensitive to the perturbations in scales k ~0.5 - O.lh Mpc -1 . 
As is seen from (2.12), estimates are very sensitive to the 
value of 2 — 7, and, so, to small scale perturbations. The first 
zero-point of autocorrelation function, r^ ^ 40/i -1 Mpc, can 
be usually found with a large uncertainty (~20 - 30 per cent) 
but its impact is reduced by the small exponent 1— 7/2 < 0.3 
in (2.12). 

For OCDM and ACDM models the impact of small scale 
matter clustering is moderate, and differences between r 5 
and T V ei are found to be ~10 per cent. The differences be- 
tween the same parameters and tt can be considered as 
a reasonable measure of simulated 'cosmic variance'. For 
these models differences between the parameter T£ calcu- 
lated for the real and redshift spaces also do not exceed ~ 
10%. For the SCDM model both and r ve i are distorted 
by the strong small scale clustering. This divergence indi- 
cates that for the SCDM model the successful application of 
methods discussed in Sec. 3 is also in question. 

The progressive growth of and as for mock catalogues 
characterizes the degree of the large scale bias between the 
spatial distribution of DM component and 'galaxies'. 



6 PROPERTIES OF WALL - LIKE 
STRUCTURE ELEMENTS 

The main basic characteristics of walls were discussed in 
DMRT for three DM and four mock catalogues mentioned 
above both in the real and redshift spaces. In this Sec. the 



wall characteristics discussed in Sec. 3 are found with the 
core-sampling technique for the same simulations and the 
same samples of walls. 



6.1 Selection of wall-like structure elements 

The sample of wall-like structure elements was selected 
with the two-parameter method described and exploited in 
DMRT. It identifies the wall-like structure elements with 
clusters found using a threshold linking length, l t hr, and 
a threshold richness, N t hr- As usual, the boundary of the 
clusters is denned by the threshold overdensity, S t hr, which 
is connected with the threshold linking length by 

3 (6.1) 



Sthr = 



(n) ^{n)l 3 thr - 



The threshold richness, Nthr, restricts the matter fraction, 
f w , associated with walls. 

The main characteristics of these samples both in real 
and redshift spaces are listed in Table 2. The values of f w m 
0.4 - 0.45 are consistent with the theoretically expected 
and observed matter fraction accumulated by walls (DD99, 
LCRS1, LCRS2). The analysis performed in DMRT shows 
that for the low density models the main characteristics of 
such wall-like elements are similar to the observed charac- 
teristics of superclusters of galaxies (Oort 1983 a,b; LCRS2; 
DURS). 



6.2 DM walls in the real space 

The analysis of DM catalogues in the real space is most in- 
teresting as in this case we can study the clear signal from 
the gravitational interaction of compressed matter and can 
reveal and characterize statistically the matter relaxation. 
Five basic characteristics of DM walls discussed in Sec. 3, 
namely, the wall thickness, h w , the dispersions of wall ve- 
locities, a v , the velocity dispersion of matter compressed 
within walls, w w , the dimensionless surface density, q w , and 
mean separation of walls, D sep , can be found with the core- 
sampling method and can be compared with those found 
in DMRT. The surface density of walls is closely connected 
with the size of proto-walls as discussed in DMRT. 

Comparison of such characteristics of matter distribu- 
tion as T V ei listed in Table 1 and t v and r m related to the 
wall properties allows us to test the influence of small scale 
matter clustering and other random factors discussed in Sec. 
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Figure 1. PDFs of DM wall velocity, N v (v w /(v w )), in real space 
for SCDM, ACDM and OCDM models. The Gaussian fits arc 
shown by solid lines. 
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Figure 2. PDFs of the DM wall surface density, N m (q w / (q w )), 
in real (solid lines) and redshift (dashed lines) spaces for SCDM, 
ACDM and OCDM models. The fits (3.1) are shown by solid lines. 



4.3, and to find the optimal ranges of core size, Dcore, and 
of threshold richness, N m in, as well as the optimal linking 
length, lu„k ■ The results listed in Table 2 are obtained with 
the linking length lunk = 5/i~ Mpc, and are averaged over 7 
core sizes, 6/i _1 Mpc< D core < 9/i _1 Mpc, and over 7 thresh- 
old richness, 10< N min <35 . 

6.2.1 Basic characteristics of DM walls 

For all models, the dispersion of wall velocities, a v , is found 
to be the best and most stable characteristic of the evolu- 
tionary stage reached. This is the direct consequence of the 
discrimination between the wall velocity and the velocity 
dispersion of particles compressed within walls. The PDFs, 
N v , plotted in Fig. 1, are well fitted to Gaussian functions 
with the measured dispersion. 

For the OCDM and ACDM models the mean dimen- 
sionless surface density of walls, (q w ), and the amplitudes, 
r m « t v w r ve i, arc found with scatters ~10 - 15% for the 
used Nmin, D CO re, and lunk- This scatter characterizes the 
moderate action of random factors discussed in Sec. 3.2 and 
the procedure of measurement. The values of l v {qw) are con- 
sistent with estimates of the size of proto walls obtained in 
DMRT. The PDFs of the surface density plotted in Fig. 2 
are consistent with that expected form (3.1). These results 
demonstrate that for lower density cosmological models the 
Zel'dovich approximation successfully describes these basic 
characteristics of rich walls. 

For the SCDM model, the results listed in Table 2 are 
more sensitive to the method of measurement and the sur- 
face density of walls is underestimated, r m < r v w T ve i. 
This difference can be mainly ascribed to the strong dis- 



ruption of walls occurring at late evolutionary stages in this 
model. Other important factors are the faster compression 
and/or expansion of walls in transversal directions and the 
existence of richer halos of evaporated particles around the 
walls mixed with infalling particles. Such a halo becomes 
richer for larger r, i.e. for the ACDM, and especially, for the 
SCDM models. 

The distribution function of wall separation, N sep , plot- 
ted in Fig. 3 is well fitted to (truncated) exponential distri- 
bution. The mean wall separation {D sep } is sensitive to the 
threshold richness N m in and to the core size D core . The sep- 
aration {D sep ) ~ 40h~ Mpc, found for the lower threshold 
richness, N min —5, and larger core sizes, D core = 9/i _1 Mpc, 
coincides with the results obtained in DMRT. It increases 
with N mirl as the number of rich walls progressively de- 
creases. For smaller D core and larger N m in some of highly 
disrupted walls are lost due to their small covering factor. 
This parameter can be found with relatively large scatter. 

Using relation (3.2) we can compare our estimates of 
(D sep ) and (q w )- For all models we have 

f dq w (0.75 - 0.9)/„ 

and the mean wall separation is probably overestimated. 

For all models under consideration, the mean wall thick- 
ness, (hw), is similar to that found in DMRT with the inertia 
tensor technique, where a wall is represented by a homoge- 
neous ellipsoid. It is about of 2 - 4 times smaller than that 
expected in the Zel'dovich approximation (3.10) what bears 
a sign of the relaxation of gravitationally bounded DM par- 
ticles within walls. 

For the OCDM model the velocity dispersion of matter 
compressed within walls is found to be similar to the mean 
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Figure 3. PDFs of DM wall separations, Nsep(Dsep/{D ae p)), in 
real (solid lines) and redshift (dashed lines) spaces for SCDM, 
ACDM and OCDM models. The exponential fits are shown by 
solid lines. 

velocity of walls and of all particles, (w w ) ~ <t v ~ cr V el- 
In contrast, for the SCDM and ACDM models the disper- 
sion (w w ) is about 30% smaller than that obtained for the 
complete walls in DMRT and the dispersions o~ ve i and a v 
discussed above. This divergence characterizes statistically 
the evaporation of high energy particles in course of the re- 
laxation of compressed matter and is reinforced by the pro- 
cedures of measurement and wall selection. The relatively 
small value of (w w ) demonstrates that in contrast to the 
clusters of galaxies the moderate degree of one dimensional 
matter compression within walls is not accompanied by an 
essential growth of velocity dispersion. 

6.2.2 Relaxation of compressed matter 

For ACDM and SCDM models the wall thickness, h w ~ (3 
- 4)/i _1 Mpc, is 2 - 3 times smaller than that expected in 
the Zel'dovich theory (3.10). So large compression of matter 
within walls means that the selected particles are strongly 
confined and, probably, relaxed. For ID matter compression 
the relaxation is expected to be weak, but in reality it is 
reinforced due to the small scale clustering and disruption 
of walls. 

The degree of relaxation reached can be characterized 
by the parameters {8} & (e) , 

Mt?}' < £ > = lfr (6 ' 3) 

listed in Table 2. Here (8) measures the mean degree of 
matter compression, and (e) is the mean kinetic energy of 
compressed particles with respect to the expectations of the 
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Figure 4. PDFs of the reduced velocity dispersion, Nu,(lo/(u))), 
for DM walls in real (solid lines) and redshift (dashed lines) spaces 
for SCDM, ACDM and OCDM models. The Gaussian fits arc also 
shown by solid and dashed lines. 

Zel'dovich theory. The function w z (t) given by (3.7) is eval- 
uated at r = t v . 

The divergence between the expectations of Zel'dovich 
theory and simulations is moderate for the OCDM model 
and becomes strong for the ACDM model as the evolution 
progresses. For the SCDM model the estimate of (8) is artifi- 
cially decreased together with {q w ). The small value of (e) ~ 
0.1 - 0.2 confirms an essential deficit of energy of compressed 
particles in comparison to that expected in the Zel'dovich 
theory. This deficit is partly enhanced by the procedure of 
wall selection, as the wall boundaries are blured, and par- 
ticles placed far from the wall center are not included into 
walls. 

In Zel'dovich theory the strong correlation of w w and 
h w with the wall richness, m w , is described by expressions 
(3.6) & (3.9). In simulations the measured linear correla- 
tion coefficients of q w , w w and h w are also ~ 0.4 - 0.5 what 
indicates that the essential mass dependence of these pa- 
rameters remains also after relaxation. To discriminate the 
regular and random variations of functions w w and h w we 
will consider the reduced wall thickness, and the reduced 
velocity dispersion, lu, which can be defined as follows: 

h w = {h w }fi Ph C, w w = {w w )fi Pw cu, (6.4) 
/i = m w /{m w ) = q w /{q w ), p h ~ Pw ~ 0.3 - 0.4, 

(C) « (w) « 1, ere « a* « 0.2. 

In all considered cases the PDFs of the reduced velocity 
dispersion within walls, N w , and of the reduced wall thick- 
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Table 2. Wall characteristics in real and redshift spaces. 
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real space 














SCDM 


2.5 


0.44 


0.74 


1.00 ±0.18 


0.39 ± 0.04 


0.58 


4.7 


3.4 ±0.4 


463 ± 18 


0.1 


0.27 


43 ± 8 


0.4 


ACDM 


1.6 


0.46 


0.85 


0.83 ±0.16 


0.35 ± 0.03 


0.39 


7.4 


4.0 ± 0.3 


387 ± 20 


0.2 


0.30 


71 ± 13 


0.4 


OCDM 


1.3 


0.40 


0.88 


0.52 ±0.06 


0.28 ± 0.02 


0.26 


2.4 


6.0 ± 0.7 


330 ± 33 


0.8 


0.44 


46 ± 8 


0.3 


mocki 


1.6 


0.43 


0.82 


0.87 ± 0.15 


0.36 ± 0.03 


0.27 


5.0 


7.0 ± 1.1 


412 ± 51 


0.9 


0.48 


84 ± 17 


0.3 


mock2 


1.3 


0.44 


0.81 


0.84 ±0.13 


0.36 ± 0.03 


0.27 


5.1 


6.6 ±0.9 


394 ± 45 


0.8 


0.47 


78 ± 13 


0.3 


mock3 


1.3 


0.45 


0.84 


0.88 ±0.11 


0.36 ± 0.02 


0.27 


5.6 


5.2 ±0.6 


354 ± 27 


0.6 


0.43 


73 ± 11 


0.4 


mock4 


1.3 


0.44 


0.86 


1.23 ±0.17 


0.43 ± 0.03 


0.28 


9.1 


4.5 ±0.4 


359 ± 20 


0.5 


0.39 


78 ± 14 


0.5 


mock4 


1.3 


0.48 


0.87 


1.35 ± 0.19 


0.45 ± 0.03 


0.27 


7.5 


5.6 ± 0.5 


379 ± 26 


0.7 


0.40 


81 ± 12 


0.5 



SCDM 


2.5 


0.41 


0.88 


0.84 ± 0.08 


ACDM 


2.1 


0.45 


0.88 


0.76 ±0.13 


OCDM 


1.3 


0.44 


0.77 


0.56 ±0.10 


mocki 


1.2 


0.43 


0.82 


0.89 ±0.16 


mock2 


1.5 


0.43 


0.83 


0.86 ±0.13 


mock3 


1.8 


0.44 


0.86 


0.85 ±0.09 


mock4 


1.8 


0.45 


0.84 


1.10 ±0.13 


mock4 


1.3 


0.46 


0.82 


1.22 ±0.13 



redshift 


space 




0.36 ±0.02 




1.8 


8.5 ±0.8 


0.34 ± 0.03 




2.6 


7.2 ±0.6 


0.29 ± 0.03 




1.4 


11.2 ± 1.3 


0.37 ± 0.03 




2.7 


13.3 ± 2.3 


0.36 ± 0.03 




2.7 


12.2 ± 1.9 


0.36 ± 0.02 




2.9 


9.5 ± 1.1 


0.41 ± 0.02 




4.1 


8.8 ±0.9 


0.43 ± 0.02 




3.8 


10.1 ±0.9 



245 ± 23 




0.75 


45 ± 9 


0.3 


207 ± 17 




0.48 


63 ± 11 


0.4 


323 ± 36 




0.58 


49 ± 8 


0.3 


385 ± 66 




0.71 


87 ± 18 


0.3 


353 ± 54 




0.70 


86 ± 25 


0.3 


276 ± 32 




0.70 


83 ± 17 


0.3 


254 ± 25 




0.65 


87 ± 21 


0.4 


291 ± 26 




0.61 


89 ± 16 


0.4 



Here 5 t ^ r is threshold parameters of clusters, f w is the fractions of all particles forming the selected walls, and f cr is the fraction of particles 
belonging to walls in both real and redshift spaces. Parameters r m ,T„ are the amplitude of perturbations as given by (3.1) & (3.4). The 
other quantities are explained in the text. Averaging was performed over 7 core sizes, 6/i _1 Mpc< D core < 9/i _1 Mpc, and over 7 threshold 
richness, 10< N m i n <35 for the linking length lunh = 5/i~ 1 Mpc. 



ness, N^, can be roughly fitted to Gaussian functions. The 
PDFs are plotted in Fig. 4 for all three models. 

These results show that due to the strong relaxation of 
compressed matter the correlations between the considered 
characteristics of walls predicted by the Zel'dovich theory 
in Eqs. (3.6), (3.8) & (3.10) are replaced by relations (6.4) 
which are also universal. 



6.3 DM walls in the redshift space 

If the analysis of wall characteristics in the real space al- 
lows to reveal the influence of gravitational interaction of 
the compressed matter, then a similar analysis performed 
in the redshift space reveals the influence of random veloci- 
ties on the observed characteristics of the large scale matter 
distribution. 

In the redshift space the analysis of wall characteristics 
was performed for samples of walls selected as described in 
Sec. 6.1 . As was shown in DMRT, in low density cosmologi- 
cal models the main characteristics of these walls are similar 
to the observed characteristics of superclusters of galaxies. 
The determination of wall characteristics and their correc- 
tions are discussed in Sec. 4.3. The wall parameters were 
found in the same ranges of D core and N min as in the real 
space for kink = 5/i _1 Mpc. 

The main results are listed in Table 2 and are plotted 
in Figs. 2-4. 

6.3.1 Walls in the real and redshift spaces 

The samples of walls selected in the real and redshift spaces 
are not identical with each other due to influence of ran- 
dom velocities of particles. This difference can be suitably 
characterized by the fraction of the same particles assigned 
to walls in both spaces. Here this fraction was defined as 



a ratio of number of the particles, N corn , to the number of 
particles assigned to the selected walls, N w . For all models 
under consideration this fraction, listed in Table 2, is 

U = N com /N w ~ 0.8 - 0.9. 

Small variations of number of particles, N w , assigned to 
walls in the real and redshift spaces lead to these variations. 

These results indicate that the influence of high ran- 
dom velocities generated by the small scale wall disruption 
and the matter relaxation moderately distorts the sample 
of walls selected in the redshift space. More strong devia- 
tions between such wall parameters as the wall thickness 
and degree of matter compression, measured in the real and 
redshift spaces, are caused by the redistribution of matter 
within walls and the procedure of measurement rather than 
by the incorrect wall identification. The impact of these fac- 
tors rapidly increases with r m and becomes extreme for the 
SCDM model. 

These deviations can be sensitive to the code used for 
simulation (see, e.g., discussion in Splinter et al. 1998). For 
example, in the P 3 M code, these variations may increase due 
to larger velocities of compressed matter generated there. 

6.3.2 Basic characteristics of DM walls 

For all three models the mean surface density of selected 
walls listed in Table 2 is similar to that found in the real 
space. This fact shows that the artificial growth of matter 
concentration within walls discussed in Sec. 3.2 is effectively 
suppressed by the influence of the velocity dispersion and 
the procedure of wall selection, and the relation (3.1), as 
before, connects the mean surface density of selected walls, 
(q w ), with the amplitude, r m . Variations of (q w ) and T m with 
D cor e and N m in are shown in Table 2 as a scatter of these 
parameters. The PDFs N m plotted in Fig. 2 are also similar 
to those found in the real space. 
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The mean wall separation is consistent with the es- 
timate found in real space and, as before, for all models 
fdq « (0.75 - 0.9)/*,. The PDFs N sep plotted in Fig. 3 are 
also similar to those found in the real space. 

In the redshift space the used method of wall identifi- 
cation selects mainly particles with a small relative velocity 
what essentially restricts the measured velocity dispersion 
within walls and the wall thickness. Results listed in Table 2 
show that only for the OCDM model, the velocity dispersion 
of compressed matter is consistent with the values found in 
the real space and in DMRT. For ACDM and SCDM mod- 
els they are even smaller than those found in the real space. 
The measured wall thickness is now linked with the velocity 
dispersion by the relation (3.16). 



6.3.3 Characteristics of matter relaxation 

In the redshift space walls are less conspicuous than in the 
real space but, even so, for all three models the mean over- 
density, (S), listed in Table 2, differs from the estimates 
based on the Zel'dovich theory (3.7). As in the real space, 
the velocity dispersion in the redshift space is strongly cor- 
related with the surface density of walls, what is described 
by the relation (6.4) with an exponent p w « 0.5. The PDFs 
of the reduced velocity dispersions, AT W , plotted in Fig. 4, 
demonstrate some excess of particles with lower w, but it can 
also be roughly fitted to a Gaussian function with (lu) wl 
and dispersion « 0.4. This dispersion is about two times 
larger than that in the real space. 

These results show that, although in the redshift space 
walls are not so conspicuous as in the real space, in the range 
of 'time' ~ 0.2 < r < ~ 0.5, the relaxation of compressed 
dark matter can be directly recognized with these methods. 



6.4 Walls in mock catalogues 

The analysis of mock catalogues characterizes how the con- 
sidered simple model of large scale bias influences the mea- 
sured wall properties. These catalogues were investigated 
also both in the real and redshift spaces. The analysis was 
performed for 10 values N min (15< N min <60) and for 
7 values of the core radius D core (7/i _1 Mpc< D core < 
10/i _1 Mpc) using a linking length lu„k = 5/i _1 Mpc. The 
main results averaged over these N m i n and D core are listed 
in Table 2. 



6.4.I 'Galaxy' walls in the real space 

In the real space for all mock catalogues the parameters 
t v , (h w ) and (e w ) are similar to those found for the basic 
OCDM model. The velocity dispersion of 'galaxies' within 
walls, (w w ) ~ a v ~ cr ve i, exceeds that found for the basic 
model by about of 20 - 30%. These variations can be at- 
tributed to the preferential identification of 'galaxies' in the 
central high density regions of walls, where the relative ve- 
locities of DM particles are also larger than the mean values. 
The wall thickness and the velocity dispersion of 'galaxies' 
can be reduced and turned into dimensionless quantities in 
the same manner as in Eq. (6.4), and the PDFs for the re- 
duced wall thickness and velocity dispersion within walls, N( 
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Figure 5. PDFs of wall surface density, N m (q w / {q w )), wall 
separation, N sep (D sep / (D sep )), reduced velocity dispersion, 
N w (w/(w)), and velocity of walls, N v (v w /a v ), for the mock^ cat- 
alogues in the real (solid lines) and redshift (dashed lines) space. 
The same fits as in Figs. 1-4 are plotted as well. 

and Nuj, are also similar to Gaussian functions. The PDFs 
7V W are shown in Fig. 5 for the mock.4 catalogue. 

As was expected, the mean surface density of walls, 
(q-w), exceeds that found for the basic OCDM model, and 
this excess progressively increases together with the biasing 
factor used. This excess can be considered as a suitable mea- 
sure of the bias. This means that to characterize this bias the 
difference between r m and t v and/or between r m and other 
amplitudes measured for the same catalogues can be used 
together with the autocorrelation function. The growth of 
q w leads to a proportional growth of 8, as the wall thickness 
is only weakly distorted. 

The large scale bias increases the contrast between 
richer and poorer walls what is seen as an essential growth 
of the mean wall separation. In all mock catalogues (D sep ) is 
about two times larger than in the basic OCDM model. The 
growth of both (q w ) and (D sep ) do not distort the relation 
between f w and fd q . 

6.4-2 'Galaxy' walls in the redshift space 

In the redshift space the fraction of the same particles as- 
signed to walls both in the real and redshift spaces becomes 
fcr ~ 80 - 85% (Table 2), what explains the similarity of 
parameters (q w ) and r m listed in Tables 2 for both cases. 
The expected growth of wall richness in the redshift space 
according to (3.13) is not found, and the surface densities of 
walls, (q w ), are, in the range of errors, the same both in the 
real and redshift spaces. This fact shows that for 'galaxies' 
the expected growth of wall richness in the redshift space 
is suppressed even more strongly than for DM component 
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due to the relaxation of compressed matter. Then Eq. (3.1) 
describes correctly the time dependence of the mean wall 
surface density. 

The parameters (h w ) & (w w ) for 'galaxy' walls are sim- 
ilar to those found for the underlying OCDM model. The 
difference of (w w ) found for the same samples of walls in 
the real and redshift spaces and a slow decrease of (w w ) for 
stronger biased models can be assigned to the loss of a small 
fraction of particles with large velocities, what demonstrates 
the sensitivity of these functions to the method of wall iden- 
tification. 

In the redshift space we have not so a reliable indepen- 
dent estimator of the amplitude as r v . There the bias is seen 
as a relation of the amplitudes r m > . This makes it diffi- 
cult to quantitatively estimate the relatively moderate large 
scale bias in observed catalogues because both r m & arc 
sensitive to the bias. 

This discussion shows that the simple algorithm used 
in DMRT for the 'galaxy' identification does not essentially 
distort the basic characteristics of simulated walls, and a 
stronger bias can be seen as an excess of the surface density 
of 'galaxies' relative to that found for the DM component 
in the basic model. At the same time the mean velocity 
dispersions of both the DM component and the 'galaxies' 
assigned to walls, (w w ), tends to be smaller than <j v & r v 
and other characteristics of the amplitude of perturbations. 



7 SUMMARY AND DISCUSSION 

In this paper we continue the investigation of large scale 
matter distribution and processes of large scale structure for- 
mation and evolution. Some aspects of these problems were 
discussed in our previous papers (LCRS-1, LCRS-2, DURS, 
DMRT, Miiller et al. 1998) where the 3-dimensional analy- 
sis of the observed and simulated large scale structure was 
performed with the core-sampling and the Minimal Span- 
ning Tree techniques. Another approach to this problem, 
based on the percolation technique, was discussed in Sahni 
et al. (1994), Shandarin & Yess (1998) and Sathyaprakash et 
al. (1998). The statistical description of structure formation 
and evolution based on the Zel'dovich theory of nonlinear 
gravitational instability can be found in Lee & Shandarin 
(1998) and DD99. 

Here we direct our attention to the physical aspects of 
the process of wall formation, what implies a more detailed 
discussion of the properties of DM walls in real space. The 
simulations described and investigated in DMRT are used 
to test the theoretical expectations, to estimate the influ- 
ence of small scale clustering and relaxation of compressed 
matter and other random factors, and to examine the power 
of the statistical methods used to describe the large scale 
matter distribution. Three cosmological models, at different 
evolutionary stages, were analyzed in the same manner, and 
the comparison of results obtained for these models allows 
us to estimate the properties of walls at various r. 

In the redshift space the influence of small scale clus- 
tering and large velocity dispersion of compressed matter 
noticeably distorts some characteristics of the walls. These 
distortions appear also in the considered mock catalogues, 
and can even be enhanced by the possible large scale bias 
between the spatial distribution of DM and galaxies. 



Some of these results may depend on the code used for 
the simulations (see, e.g., the discussion in Splinter et al. 
1998), and they should be checked with simulations employ- 
ing a code with higher spatial resolution. 

7.1 Identification of walls 

The core-sampling approach described in Sec. 5 allows us to 
characterize, in more details, the matter distribution along 
the sampling core and to estimate the uncertainty in mea- 
sured properties of wall-like condensations introduced by the 
influence of velocity dispersion and small scale clustering. 
The influence of these random factors is demonstrated by 
comparing results obtained with various D core , N min & lu„k- 

Results presented in Sec. 6 show that some fraction of 
the early compressed matter has subsequently evaporated 
due to relaxation processes. These DM particles together 
with the infalling matter form an extended halo around the 
walls and, therefore, it is difficult to separate the walls from 
the background. The same problem is met by the correct 
definition of boundaries of galaxies and clusters of galaxies. 
It was also discussed in the DMRT, LCRS-2 and DURS, 
where the methods of wall selection, described in Sec. 6, were 
applied to simulated DM and observed galaxy distributions. 

The central high density part of walls is reliably selected 
in all the cases, but various definitions of the wall boundaries 
can noticeably change the measured characteristics of walls. 
To provide more objective comparisons of wall characteris- 
tics the same dimensionless parameters f w and 8 t hr should 
be used for identification of walls in different catalogues and 
simulations. 

7.2 DM walls in the real space 

7.2.1 Measured characteristics of walls 

The results presented in Sec. 6 show that the core-sampling 
approach can be successfully used for the investigation and 
description of the large scale matter distribution and the 
wall-like matter condensations. It allows to estimate the sur- 
face density, thickness, velocity dispersion and other basic 
parameters of DM walls corrected for the influence of ran- 
dom curvature and shape of walls. These parameters differ 
from those obtained in 3D space with the Minimal Span- 
ning Tree and inertia tensor methods, and these methods 
suitably complement each other. 

The measured wall characteristics can be compared 
with predictions of the Zel'dovich theory what reveals the 
influence of relaxation of compressed matter on the proper- 
ties of walls and allows to correct the theoretical expecta- 
tions. The small scale clustering of compressed matter and 
the wall disruption lead to noticeable variations of measured 
wall characteristics for different parameters of the sampling 
core. These variations are not so large for the low density 
models, but they increase rapidly with r. 

The dimensionless surface density of walls, q w , is closely 
connected with the size of proto- walls as discussed in DMRT, 
LCRS2 and DURS. The high surface density of walls, q w > 
0.6, found above even for the low density models, demon- 
strates that processes of strong nonlinear matter evolution 
occur at a typical scale of ~ q w l v ~ (15 - 25)ftT 1 Mpc. This 
evolution is correctly described by the Zel'dovich theory. 
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This characteristic is sensitive to the basic cosmological pa- 
rameters, Q m & h, what allows us to select the class of more 
perspective models for further investigation. 

7.2.2 Relaxation of compressed matter 

The problem of relaxation of compressed matter is now in 
the forefront, and the obtained results allow to begin discus- 
sion of the statistical characteristics of this relaxation. The 
analysis performed in the real space is more important for 
the discussion of the basic physical processes which occurred 
during the formation of wall-like matter condensations, such 
as the small scale matter clustering and the relaxation of the 
compressed matter. These processes generate the large ve- 
locity dispersion within walls and lead to the evaporation 
of high velocity particles. Thus, in all these cases an es- 
sential deficit of energy in DM walls as compared with the 
expectations of the Zel'dovich theory - ~ (50 - 80)% and 
more - was found. The growth of this deficit with r from 
the OCDM to SCDM models demonstrates that the DM re- 
laxation becomes more and more important for later evolu- 
tionary stages, and its influence on the observed parameters 
of the large scale matter distribution becomes crucial for 
r > 0.5. 

The relaxation is seen in rich superclusters of galaxies 
such as the Perseus-Pisces (Saslaw & Haque- Copilah 1998). 
It is essentially accelerated and amplified by the small scale 
clustering of compressed matter. This clustering is clearly 
seen in observations as, for example, a strongly inhomoge- 
neous galaxy distribution within the Great Wall (Ramella et 
al. 1992). The clusters of galaxies situated within wall-like 
superclusters similar to the Great Wall and the Perseus- 
Pisces can be considered as extreme examples of this pro- 
cess. 

The merging of earlier formed structure elements is very 
important for the formation of large walls (DD99). This 
means that actually the relaxation occurs step by step dur- 
ing all the evolutionary history beginning with the forma- 
tion of first low mass, high density pancakes which later are 
successively integrated and merged to larger structure ele- 
ments. This means also that the finally reached degree of re- 
laxation and the properties of compressed matter depend on 
the (unknown) evolutionary history of the considered walls 
and, therefore, can be characterized only statistically. 

The relaxation of compressed matter destroys the tight 
correlation between the surface density and velocity disper- 
sion predicted by the Zel'dovich theory (3.6), but it gen- 
erates other correlations between the same characteristics 
described by the relations (6.4). This fact indicates that the 
properties of compressed matter are sufficiently general, and 
these characteristics can be used to improve the methods of 
wall selection and the description of wall properties. 

The velocity dispersion within walls increases gradually 
with r from the OCDM to the SCDM model. As was dis- 
cussed in Sec. 7, the particles with high velocity are gravita- 
tionally confined and occupy preferentially the high density 
central regions of walls. This fact confirms that, probably, 
these particles are relaxed and have a (quasi) stationary dis- 
tribution. This distribution is not so stationary as, for ex- 
ample, in clusters of galaxies, and it is slowly evolving due to 
the large scale matter flow along the walls and the persisting 
merging of neighboring structure elements, but presumably, 



this evolution does not significantly distort the formed mat- 
ter distribution. 



7.3 DM walls in the redshift space 

The matter condensation seen in the redshift space can be 
partly artificially enhanced by the influence of streaming ve- 
locities. The possible influence of this effect was widely dis- 
cussed over the past decade (see, e.g., Kaiser 1987, McGill 
1990 a; Davis, Miller & White 1997; Hamilton 1998; Hui, 
Kofman & Shandarin 1999) and, as applied to properties of 
absorption lines in the spectra of high redshift quasars, by 
McGill (1990b) and more recently by Levshakov & Kegel 
(1996, 1997). These tendencies are also clearly seen from 
the direct application of the Zel'dovich approximation to 
the wall formation in the redshift space as was discussed in 
Sec. 3.2. Of course, it is impossible to decide which particles 
belong to walls, but we can estimate statistically the prop- 
erties of DM walls identified in the redshift space. However, 
the influence of this uncertainty cannot be separated from 
the influence of the relaxation and other factors discussed 
above. 

For all models the comparison of DM walls selected in 
the real and redshift spaces demonstrates, that they are com- 
posed mainly from the same particles - this fraction is about 
for ~ (70 - 80)% (Table 2). This means that in both cases 
we find the same walls, and the fraction of randomly added 
or lost particles is indeed small. In spite of this, some prop- 
erties of walls in the redshift space are quite sensitive to the 
velocity dispersion and to the methods of wall identification. 
Thus, the strong growth of wall thickness - about a factor 
of 2 - confirms results obtained by Melott et al. (1998). This 
effect is quite similar to the well known 'finger of God' effect 
observed in clusters of galaxies. 

The wall surface density, q w , is most interesting, as it is 
directly connected with the basic cosmological parameters, 
Q m & h. Our analysis shows that for low density models 
- ACDM and OCDM - the measured value of q w is sim- 
ilar both in the real and redshift spaces. This means that 
the growth of the matter condensation within walls due to 
streaming velocities as predicted by the Zel'dovich theory 
is strongly suppressed by the influence of the matter relax- 
ation and the transformation of a continuous matter infall 
to a discontinuous one. Actually similar relations connect 
the fundamental wall characteristics such as q w and r m . 

The velocity dispersion within walls selected in the 3D 
redshift space can be noticeably underestimated what is a 
direct consequence of the method of wall selection. In the 
redshift space, particles with large velocities are artificially 
shifted to the periphery of selected walls and, so, can be 
omitted from the analysis. 

7.4 Walls in mock catalogues 

For the considered mock catalogues the influence of veloc- 
ity dispersion is enhanced by the methods used for 'galaxy' 
selection. The large scale bias increases the 'galaxy' concen- 
tration within walls and, so, increases the density gradient 
near the wall boundary. When the 'galaxies' are identified 
preferentially in the high density central parts of the walls 
(in the real space), than their velocity dispersion exceeds 
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that for the DM particles, and this excess may be as large 
as ~ (20 - 30)%. In the redshift space, the parameters of 
'galaxy' walls such as (h w ) and (w w ) are similar to those in 
the underlying DM distribution. 

The bias is clearly seen both in the real and redshift 
spaces as an excess of the mean surface density of walls. 
The comparison of parameters q w & r m found for observed 
wall-like galaxy condensations with possible independent es- 
timates of the same parameters gives us a chance to obtain 
a reasonable observational estimates of the large scale bias. 



7.5 The amplitude of large scale perturbations 

These results demonstrate again that all characteristics of 
the amplitude and evolutionary stage of large scale structure 
considered in Sees. 2 & 3 are similar, but not identical to 
each other, as they are sensitive to different properties of 
perturbations. The best and most stable measure, t v , comes 
from measurements of the velocity of structure elements. 
It is insensitive to the nonlinear evolution of perturbations, 
large scale bias and small scale clustering or relaxation of 
the compressed matter. 

The comparison of other estimates for the same param- 
eter r, namely, T ve i, tj, & r m obtained in the same simula- 
tions demonstrates their sensitivity to various natural and 
artificial factors. For the low density models - ACDM and 
OCDM - the parameters t v and r m are usually sufficiently 
close to each other, what is a direct consequence of the close 
connection of the process of wall formation with the large 
scale perturbations. The parameter r m is sensitive to a pos- 
sible large scale bias, but to reveal this factor, we need to 
have independent unbiased estimates of the same amplitude. 

The most interesting independent estimate of the am- 
plitude is T£ which is however more sensitive to small scale 
matter clustering. Thus, for the SCDM model where this 
clustering is stronger it significantly overestimates the large 
scale amplitude. It is less sensitive to the large scale bias 
than T m ■ 

Independent estimates of the large scale amplitude 
come from measurements of the CMB anisotropy. The 
COBE data are consistent with other available estimates 
of cosmological parameters and of the large scale ampli- 
tude, and therefore, tt can be considered as the best es- 
timate of the combination (2.9) of F and the amplitude. It 
can be connected with estimates of cosmological parame- 
ters fi m « 0.3, Qa ~ 0.7 obtained from observations of 
high- redshift supernovae (Perlmutter et al. 1998). Nonethe- 
less, tt should be corrected for a possible contribution of 
gravitational waves. 

The investigation of the space density of clusters of 
galaxies and its redshift evolution (see, e.g., Bahcall & Fan 
1998; Eke et al. 1998; Wang & Steinhardt 1998) seems also 
to be promising and can give the required independent mea- 
sure of the large scale amplitude. The formation and evolu- 
tion of galaxy clusters is caused by large scale perturbations, 
and their characteristics can be connected with these pertur- 
bations. But they are sensitive to the thermal evolution of 
clusters and, moreover, are related to only ~ (10 - 15)% of 
matter accumulated by the clusters. This means that they 
are not free from random variations what is seen, in par- 
ticular, as the well known variations of the autocorrelation 
function with the cluster sample. 



The critical discussion of available measurements of cos- 
mological parameters (Wang et al. 1999; Efstathiou 1999) 
shows that in spite of a large progress reached during last 
years, we do not have yet a reliable unbiased estimate of 
these parameters, and these data should be tested with re- 
spect to possible random large scale variations. The appli- 
cation of the discussed methods to large observed redshift 
surveys can help to achieve this goal. 
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Appendix A 

Dynamical characteristics of walls 
in the Zel'dovich theory 

The results obtained in DD99 allow us to discuss in more 
details dynamical characteristics of walls predicted in the 
Zel'dovich theory. The comparison of these expected and 
actually simulated characteristics reveals the influence of in- 
teraction and relaxation of compressed matter. 

Following DD99 we define the wall formation as the 
intersection of two DM particles with different Lagrangian 
coordinates, qi and q 2 . The difference of these coordinates 
measures the size of the pancake. Using the basic relations of 
the Zel'dovich theory (2.2) and (2.3), linking the Lagrangian 
and Eulerian coordinates and velocities of particles, we ob- 
tain the coordinate and velocity of a wall as a whole (DD99): 



l v H(z) 
1 + z 



qi q2 
|qi - q2| 



q c -r(z)(l + (3) 



A<f>(q u 



Qc = 



|qi + q2 



21 1< 



q w = 



qi - q2 



l [! 



where A&(q w ) is the random difference of the dimensionless 
gravitational potential over the wall. It is convenient to in- 
troduce the relative normalized Lagrangian coordinate of a 
particle within a wall, 



q c + 0.5q u 



Using the coordinate ■& we will describe the relative position 
and velocity of the infalling particle with the Lagrangian 
coordinate q p or by the functions: 
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v in f = nv - v w = —u{z)Q.^q w d + H(z)(l + (3)r inf , (A. 2) 

u(z) = H(z)l v (3(z)(l + z)-\ 

Here S = nS is thee random dimensionless longitudinal dis- 
placement of a particle from its unperturbed Lagrangian 
position introduced by (2.2). 

For Gaussian initial perturbations the PDF of the ran- 
dom function r in f is also Gaussian, and the mean value 
and dispersion of r in f should be found using the conditional 
characteristics of functions S and A$ taking into account 
that a wall is formed in the point r — r w (DD99). In this 
case for walls with (}„ <1 we have: 
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and (?"?„*) is independent from d. This means that both 
random functions, 

rt„f & v in f + u(z)0.5q w d = H(z)(l + /3)r in f 

are also independent from Hence, for the thickness, h w , 
of a wall with the surface density q w , and for the velocity 
dispersion within such a wall we have 
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Here the wall thickness is normalized by the thickness of a 
homogeneous slice. 
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